Stability of quantum motion in regular systems: a uniform semiclassical approach 
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We study the stability of quantum motion of classically regular systems in presence of small 
perturbations. On the base of a uniform semiclassical theory we derive the fidelity decay which 
displays a quite complex behaviour, from Gaussian to power law decay t~" with 1 < a < 2. 
Semiclassical estimates are given for the time scales separating the different decaying regions and 
numerical results are presented which confirm our theoretical predictions. 
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Stable manipulation of quantum states is of impor- 
tance in many research fields such as in quantum in- 
formation processing and in Bose-Einstein condensation. 
A measure of the stability of quantum motion is the 
so-called fidelity or quantum Loschmidt echo 1], which 
characterizes the stability of quantum dynamics under 
small perturbations of the Hamiltonian that may derive 
from static imperfections or from interaction with an ex- 
ternal environment. From a more general point of view 
it is interesting to understand the behavior of fidehty in 
relation to the dynamical properties of the system. 

While fidelity decay in classically chaotic systems has 
been extensively studied 0, i, i, SM S B i, [iSlllL the 
situation in regular systems is much less clearjsl Hol [iTI . 
m, El, Q m, Ili, [13, [H, m and only for the particular 
case of vanishing time-average perturbations a clear un- 
derstanding has been achieved [I3 • Moreover, for a gen- 
eral perturbation and for a single initial Gaussian wave 
packet, a Gaussian decay has been predicted in some time 
interval which is not exactly specifiedQ. On the other 
hand, numerical investigations show a much more rich 
behavior of the fidelity decay ranging from power law to 
exponential, up to a Gaussian decay, depending on initial 
conditions, perturbation strength, and time interval^l5|. 
In addition, a somehow unexpected regime in which the 
fidelity decay in regular systems is faster than in classi- 
cally chaotic systems has been found in Q . All the above 
calls for a theory which can account for these diverse an- 
alytical and numerical findings. 

We would like to draw the reader's attention to the 
fact that in the general theory of dynamical systems, in- 
tegrability is the exception rather than the rule and it is 
therefore extremely rare. However, the theory we present 
in this paper applies also to integrable regimes of systems 
with divided phase space in which both chaotic and in- 
tegrable components are present and this is indeed the 
typical situation. Moreover as shown in [20| . for a proper 
operability of a quantum computer it is desirable to re- 



main below the border for transition to quantum chaos, a 
situation which is likely to correspond to quasi integrable 
behavior. Finally, some quantum algorithm can have a 
phase space representation and, in particular, the Grover 
algorithm can be interpreted as a simple quantum map 
which turns out to be a regular map 21]. 

The above considerations motivate our interest in the 
stability of integrable motion. In this letter, we develop 
a uniform semiclassical approach to the fidelity decay in 
regular systems and we provide a unified description to- 
gether with the corresponding time scales. Numerical 
computations confirm our analytical estimates. 

Quantitatively, the fidelity for an initial state |5'o) is 
defined as M{t) = |rn(t)p, where 

m{t) = (*o|exp(ii/t/?i)exp(-ii/ot/?i)|*o)- (1) 

Here Hq and H = Hq + eV are the unperturbed and 
perturbed Hamiltonians, with e a small quantity and V 
a generic perturbing potential. 

Consider an initial Gaussian wave packet in a 2d- 
dimensional phase space, centered at (ro,po), 

V'o(ro) = K')"'^'exp [zpo • ro/n - (ro - ?o) 7(2^')] • 

(2) 

For a sufficiently narrow, initial Gaussian packet, the 
semiclassical approximation to the fidelity amplitude is 
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dpo exp 



,(3) 



where Wp = h/^ and AS* is the action difference between 
the two nearby trajectories of the two systems H and Hq 
starting at (po,?o)[X]- We mention that, for not very nar- 
row initial Gaussian packets, Eq. ^ may still hold with 
a redefinition of Wp'sil. For more general initial states, 
the fidelity amplitude can be expressed in terms of the 
Wigner function of the initial state in a uniform semi- 
classical approach (To| . 
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FIG. 1: Comparison between the numerically computed fi- 
delity decay ( black circles and open squares) and the semi- 
classical prediction Msci(t) in Eq. ([§]), in the kicked rotator 
model with k = 0.3,(7 = 1.5, iV = 2", and = n/2Q. Cen- 
ters of the initial Gaussian packets are: (ro,po) = (l-Svr, O.Gtt) 
for the circles, and (1.27r, 0.27r) for the squares. The solid and 
dashed curves represent the semiclassical predictions, respec- 
tively, with kp evaluated numerically in the classical systems. 
Inset: AS/2tt vs po/2n, for ro = 1.2n at t — 20. The aver- 
age slope of AS in the monotonically increasing part is much 
larger than that in the central part, therefore, the fidelity 
decay at po = 0.27r is faster than that at po — O.Btt. The os- 
cillations of AS in the central region imply larger second and 
higher order terms in the expansion of AS around po = O.Gtt, 
hence, larger deviations of the fidelity from Msci{t). 



The action difference can be calculated in the first or- 
der classical perturbation theory: 



AS": 



dt'V[r{t'),p{t')] 



(4) 



with V evaluated along one of the two trajectories. Equa- 
tions ^ and dl]) give quite accurate predictions even 
for relatively long times, much more accurate than that 
usually expected for a first order perturbation treatment 
0, The reason for the unexpected accuracy is ex- 
plained in by making use of the shadowing theorem, 
the two trajectories for AS" being in fact the so-called 
shadowing trajectories of the two systems with slightly 
different initial conditions. 

Equation ^ shows that the behavior of fidelity is 
mainly determined by two factors: (i) AS* as a function 
of Po and t; (ii) the Gaussian term, which specifies an 
effective window in the po space with size Wp. For sim- 
plicity, in what follows, we consider kicked systems with 
d — 1 and set the domains of r and p to be [0, 27r). 

For an integrable system or a regular region of a mixed 
system, at least locally there exist action-angle variables 
(0,1), connected to the variables {r,p) by a canonical 



FIG. 2: The time scale n versus po, for the same parameters 
as in Fig. [T] The circles give the values of n calculated by 
the first time at which \[M{t) - Msci{t)]/M{t)\ > 0.1 [see 
Eq. (|8])]. The solid curve is the semiclassical estimate given 
by aWp\kpp(Ti)\ = 1, with kpp{t) numerically computed in the 
corresponding classical system. 



transformation. The integrand in Eq. Q can then be 
written as 



Vt = V[I{t),e{t)], with e{t) ^00 + lyt, 



(5) 



I{t) ~ /o, and J/ = dH{I)/dI. 

The main features of AS" as a function of t and pq can 
be seen by substituting Eq. ^ in Eq. ^ , replacing by 
(j) = vt, and noting the periodicity of the angle 6. This 
gives 

1 /-^^ 

AS ~e[Uit + Sf), where Ui = — V{I,9)de, {6) 

27r Jo 



Sf 



V{1, 00 + (j))d(j) - bUi 



(7) 



Here h = vt — 27rnt, where nt is the integer part of vijl-K. 
Therefore, for a fixed poj ^ is a sawtooth- type function 
of time oscillating between and 27r with a frequency 
v/2ti. As a result Sf oscillates correspondingly, hence 
AS* oscillates around its linearly increasing part elJit. 
On the other hand, for a fixed t, Uj as a function of po 
changes almost linearly in the neighborhood of po while 
Sf oscillates with a frequency v't/2Tr, where v' = dv/dpo. 
Therefore the average slope of AS" with respect to po is 
eU'jt, where Uj' = dUj/dpo. (We use a tilde above a 
quantity to indicate its value taken at the center (ro,po)-) 
We first study the fidelity decay for times t < ti, where 
Ti is the time scale before which the right hand side of 
Eq. ^ can be calculated by a linear approximation to 
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FIG. 3: Solid curve: Fidelity decay for a = 0.2, N = 2^^ , 
k = 0.3, = h/20 and (ro,po) = (1.27r, O.Stt). Dashed curve: 
the prediction of Eq. (HH) with c = 1 and f7| ^ -0.019, 
U'/ ~ 0.042 numerically computed from the classical system. 
Dotted curve: the Gaussian decay c~^"'p''^^''^^^. (For smah 
values of a, the difference between this decay and Mgci in 
Eq. is small before n.) The dashed-dotted line shows the 
approximate, intermediate, exponential decay. Inset: Long 
time decay for the same case with ro — O.Gtt (every 500 steps 
shown). The dashed line gives the decay. 



AS* with respect to pq — pq. This gives 
1 



Msci(0 ^ exp 



'-{(TWpkpf 



with kp = — - — . (8) 
e dpo 



where a = e/fi is the strength of perturbation. The ex- 
plicit dependence of kp on time t can be calculated by 
using Eq. The leading terms give 



kp ~ {U'j + C/e)i, where Ug 



t dpQ 



(Ft -[//)-. (9) 
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Since Uq oscillates in time around zero, the fidelity has, 
on average, an initial Gaussian decay with a rate de- 
pending on initial conditions, with larger values of \U'j\ 
implying faster decay. 

To test the above predictions, we consider the kicked 
rotator model, H = \-<j? -I- fco cosr b(t — nT). The 

quantized system has a finite Hilbert space with dimen- 
sion N. We take fco = fciV/27r and T = 2n/N, with k = 
fcoT independent on N. The classical limit corresponds to 
N oo. The one period quantum evolution is given by 
the Floquet operator U = exp[— exp[— zfco cos(f)], 
and is numerically computed by the method of fast 
Fourier transform. Here 2'k/N serves as an effective 
Planck constant. The inset of Fig. [T] shows an exam- 
ple of AS* vs Poj in which the values of \U'j\ are large at 
the borders while are quite small in the central, oscillat- 
ing region of p^. In agreement with our theory, under the 



FIG. 4: Fidelity decay for several values of a with k — 0.3, 
iV = 2^^ ^2 = ;!./20, and (Fo,Po) = (1.27r, 0.287r). For t 
around 100, the rate of fidelity decay increases for < cr < 
0.2, then decreases for 0.2 < ct < 0.4, and then increases 
again. This is in agreement with the semiclassical prediction 
that the decay rate /3 oscillates with a period li/' /U'j\ « 0.4 
(see text). The values U'j ~ 1.85 and i?' ~ 0.76 are obtained 
from the classical motion. Notice that the transition time 
from Gaussian to power law decay of fidelity is proportional 
to 1/(7 (see text). Here the power law decay is visible for 
a = 0.4 and 0.5 (the dashed straight line gives the l/t^"^ 
decay) . 



perturbation fc — s- fc -|- e, the fidelity for initial states ly- 
ing in the two regions of po has a quite different decaying 
rate, as seen in Fig. [1] The agreement with the theory 
is particularly good in the case where the oscillations of 
AS* arc not too strong. 

The time scale ti can be estimated by the time at 
which the second-order term in the Taylor expansion of 
AS/h is order of one at the point po = pq + Wp , i.e. 
aWp I kpp (t) I ~ 1 for t — ti, where 



kpp (t) 



dkp 
dpo 



i^'Y dV , 
V d9 



0{t). 



(10) 



A numerical check of this prediction is presented in Fig.[2l 

When the term dominates, ti cx ^Xjaw^. 

Beyond ti, a direct analytical computation of the fi- 
delity is difficult since higher and higher order terms in 
the Taylor expansion of AS* with respect to po—po need to 
be considered. We take therefore the following approach: 
we divide the domain of po into segments separated by 
points poj (with poo = Po), in such a way that Sf com- 
pletes one oscillation period within each segment. (We 
will use the subscript j to indicate quantities taken at the 
point Poj). It is easy to see that {poj+i — Poj)Vjt ~ 2tt 
and therefore the number of segments increases linearly 
with time t. The quantity msc{t) in Eq. ^ can now 
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be written as a sum of contributions rrijlt) of different 
segments, namely, niscit) = J^j nijit). 

Let us introduce the time scale r^, at which Sf com- 
pletes one oscillation within the window Wp = {po — 
Wp,Po + Wp). From Eq. (O, « Tr/{i''wp). For t > r^, 
there are many segments within Wp hence, within each 
segment, the variation of the Gaussian term on the right 
hand side of Eq. ([3]) is negligible. Therefore one may 
write 



mj{t) 



-(poj-po)^/wl 
^/nwpv'jt 



„*AS,/n^(^)^ (11) 



Fit) 



f ■ i-0(t)+4> _ 
d(t)exp{ — V{I,0')d9' 
I Je(t) 



(12) 



Indeed for poj not far from pQ, Fit) is independent on 
j and since |i^(i)| does not decay with time we will not 
consider it further. 

First we notice that equation (|lip predicts a plateau in 
the fidelity decay, when the change of ASj/h within the 
window Wp is negligible. The plateau disappears when 
aWpU'jt TT. If, for example, Uj = for the system Hq, 
then Uj oc e for the system _ff , and the plateau will end at 
a time proportional to h^^^e^'^ for ^ cx h^^^ in agreement 
with the result of 12,]. 

Due to the Gaussian term on the right hand side of 
Eq. (|lip . the main contribution to vriscit) comes from 
PQj not far from p^. For these pQj, v'^ is almost a con- 
stant, as well as the value of 5'/, hence, apart from a 
common phase, AS'j in Eq. (jlip can be approximated by 
ellijt. For a sufficiently large number of segments, the 
sum '^jTTijit) can be replaced by an integral over pq. 
Then, setting q = pq — po, we have 



ruscit) 



Fit) 



dqe 



~1 /w-^iiytUi 



2'k^/t:Wp J _ 

Expanding Ui to the second order terms in we have 



(13) 



M,e2(0 



2c 



: exp 



-f- iw^aU'/t)^ [4 + iwlaU'/t)^_ 



-2iwpaU[tf 



(14) 
1 for 



where Uj = dU'j/dpo and c a constant with c 
sufficiently small a. 

From equation ([T^ it is seen that for iw'^alljt)'^ <?C 4, 
the fidelity has a Gaussian decay g^-i^v'^Vjt)- /2 ^ which 
agrees with the one given in for weak perturbations. 
On the other hand, if iwpalljt)'^ ^ 4, i.e. for larger time 
t, Msciit) in Eq. has a power law decay 1/t. In the 
transition region between the two decays, the fidelity may 
have an approximate exponential decay (see Fig. [3]) which 
can explain the exponential-like decay found numerically 
in [ip. 

With further increasing time t, higher order terms in 
the Taylor expansion of Uj will become important and 
will modify the 1/t decay. In order to evaluate the ef- 
fect of higher order terms, we divide the interval Wp 



into subintervals labelled by /, in such a way that inside 
each of them the linear approximation of Ui can be used. 
Thus their width must be of the order 5q ~ 1/ yJ\Uj\(jt 
and their number in the region Wp is Im oc ^/ai. For 
sufficiently long time t, Im ^ 1 and the width of each 
subinterval is much smaller than Wp. As a result, the 
Gaussian term in Eq. (jlSp can be regarded as constant 
within each subinterval. Then, the integral (jl3p reduces 
to niscit) ~ ^iniiit) — c/e*'^''-*^'''^; with some co- 
efficients ci and phases 4>iit). The detailed behavior 
of this sum depends on the properties of the function 
Uj. However, one can give estimates in some limiting 
cases: (i) Since | X)/ Qe^"^'^*^ I < EJc/| oc iat)^^^, the 
slowest decay is 1/^/ai. (ii) The fastest decay is ob- 
tained when I Qc*'^'^*-'! is a constant (which may hap- 
pen due to mutual cancellation of phases (j^i). In this 
case, nisc ''-^ 1/at. (iii) In the case of random phases (pi, 
one has | c/e^'^'^*) | oc (cr^)^^^ then, nisc ~_l/(c^^)3/^ 
which coincides with the result given in Ref. [l3| for av- 
eraged fidelity. (For an analysis of Gaussian and power 
law decay of fidelity, which is based on statistics of action 
difference, see Ref. ^10i]c). 

Therefore, in general, the decay of M(t) has the power 
law dependence iet/h)~°' with 1 < a < 2. This is in 
agreement with numerical results in [Tsj . as well as with 
our extensive numerical simulations (see, e.g., Fig. 4 and 
the inset of Fig. 3). The cases of a ~ 1 or 2 have been 
found quite rare in our simulations. 

It is important to remark that, contrary to chaotic sys- 
tems for which the decay rate depends on the strength 
but not on the shape of the perturbation, integrable sys- 
tems lack of such generic behavior. This is typical in 
the general theory of dynamical systems and it is due to 
the peculiarity of integrability. In the present case the 
numerical value of a depends on the particular shape of 
Ui. 

Our approach can also explain an interesting feature 
of fidelity decay in regular systems observed numerically 
in J^i, that is the fact that, in some case, the decay rate 
may decrease with increasing the perturbation strength 
a. First we note that Eq. is valid only when | AS'j+i — 
ASj\/h is small compared with tt, since it is obtained 
by replacing the sum by an integral over pq. 

Since Uj is given by (ASj+i — ASj)/ ipoj+i —poj)e^t, then 
lAS'j+i — ASjl/U ~ 2na\Ujj/i^j\, which may exceed tt for 
sufficiently large a. In this case, we can replace ASj/h 
by tpj = ASj/h — 2'Knij, where m is an integer such that 
— < TT. Then using the fact that, for pqj close 
to Pq, ipoj — pQjv't ~ 211] we find that the term cr[/j 
in Eq. (fHl) can be replaced hy (3 = \(jU'j — mv'\. With 
increasing cr, the value of (3, which gives the decay rate, 
oscillates between and \v'\/2, with a period \v' /U'j\. 
The results of Fig. 2] nicely confirm the above analysis. 

Finally a word of comment on the comparison of fi- 
delity decay in classically regular and chaotic systems. 
As shown in , in sufficiently weak perturbation regime, 
cr < (Tc oc Vh, the Gaussian decay ^-i.^p'^u'^t)'^ /2 ^j^g 
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regular case can be faster than the fermi-golden-rule de- 
cay in the chaotic case. This is not surprising as it may 
look at first sight, since in this regime, Wp<Tc oc h can 
be quite small and the fidelity may remain close to 1 for 
a time comparable to the Heisenberg time ^ f/7i, which 
can be quite long. Moreover, the above Gaussian decay in 
the regular case is followed by a power law decay, which 
is slower than the exponential decay in the chaotic case. 
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